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The  goal  of  this  pilot  project  is  to  demonstrate  unambiguous  synchronization  between  two 
nanoscale  mechanical  oscillators.  In  the  system  we  studied,  two  mechanical  resonators  are 
remotely  placed  and  linked  in  a  closed-loop  photonic  bus  waveguide,  which  can  be  treated  as  a 
photonic  cavity.  A  programmable  optical  field  is  used  to  achieve  controlled  coupling  between  the 
two  resonators  via  delayed  optical  force.  The  optical  backaction  force  induces  both  displacement 
coupling  and  velocity  coupling,  and  gives  rise  to  much  richer  dynamics  than  previous  studied 
system  where  the  coupling  occurs  in  only  one  degree  of  freedom. 

1.  Introduction  of  the  proposed  scope  of  research 

Synchronization  is  a  ubiquitous  phenomenon  where  the  phase  difference  between  free-running 
oscillators  remains  constant  due  to  the  mutual  coupling.  Besides  its  well-accepted  importance  in 
biological  sciences,  today  synchronization  is  becoming  a  powerful  tool  for  many  engineered 
systems1.  For  instance,  synchronization  is  desirable  in  situations  where  high  oscillating  power, 
strong  coherence,  or  low  phase  noise  are  needed,  such  as  lasers2,  phase-locked  loops3,  Josephson 
junction  arrays4,  and  spin-torque  resonators5.  Synchronization  also  promises  to  improve  the 
accuracy  of  time-keeping  devices6.  Since  the  observations  of  synchronization  in  pendulums7  this 
concept  has  found  its  bearings  in  science  and  engineering  due  to  its  potential  applications  in 
generating  low-noise  stable  oscillating  signals.  Nanomechanical  oscillators,  on  the  other  hand, 
are  very  appealing  as  they  simultaneously  offer  high  quality  factor  resonances,  excellent 
scalability  '  and  are  ideal  systems  for  synchronization  studies  due  to  their  highly  engineerable 
nonlinearities10. 

However,  achieving  reproducible  and  strong  coupling  in  nanomechanical  devices  remains 
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difficult  due  to  the  unavoidable  device  non-uniformity  and  weak  mutual  coupling.  This  can  be 
circumvented  by  coupling  nanomechanical  resonators  to  an  optical  cavity11.  Recently, 
synchronization  between  two  closely  spaced  micromechanical  resonators  was  demonstrated 
using  a  hybrid  optical  mode  of  two  coupled  disk  resonators  and  the  synchronization  phase-space 
was  predicted  using  the  Kuramoto  model  .  Here  we  experimentally  demonstrate  the  first 
synchronization  of  two  spatially  separated  nanoscale  radio-frequency  oscillators  integrated  inside 
an  optical  racetrack  cavity.  We  show  that  this  leads  to  a  limit  cycle  in  the  reduced  three- 
dimensional  mechanical  phase  space  (The  two  mechanical  resonators’  amplitudes  and  their 
phase  difference  )  and  that  the  dynamics  of  two  mechanical  modes  coupled  via  a  common 
optical  mode  cannot  be  captured  by  the  standard  Kuramoto  model14:  as  a  result  of  the  additional 
degrees  of  freedom  of  the  coupled  system,  slow  dynamics  appear  on  top  of  the  limit  cycle,  and 
sidebands  emerge.  These  sidebands  are  true  signatures  of  synchronized  motion  in  the  mechanical 
domain  and  are  not  to  be  confused  with  simple  nonlinear  intermodulation  oscillatory  modes. 
Their  presence  is  important  for  the  phase  noise  performance  of  synchronized  optomechanical 
oscillators,  and  could  counteract  the  common  perception  that  synchronized  states  should  always 
improve  phase  noise  performance. 

2.  Our  experimental  research  and  main  findings 

We  investigate  the  interaction  between  two  MHz-frequency  nanomechanical  resonators  that  are 
linked  in  an  optical  racetrack  (Fig.  1(a));  The  resonators  are  mechanically  isolated,  due  to  their 
large  separation  (~  80  pm),  ensuring  that  any  coupling  between  them  is  through  the  optical  field. 
The  fabrication  of  these  integrated  photonic  devices  is  readily  scalable  [1,  15],  making  this  an 
ideal  platform  for  synchronization  studies  [11,  14].  The  silicon  beams  are  slightly  buckled  and 
they  may  end  up  in  the  up  or  down  state  (Fig.  1(a)). 

The  measurement  setup  shown  in  Fig.  Id  consists  of  a  strong  pump  laser  to  create  cavity 
backaction  and  a  weak  probe  to  detect  the  motion  [13].  First  the  two  nanomechanical  resonators 
are  characterized  with  the  pump  laser  off,  so  that  the  backaction  is  small.  When  the  two 
resonators  are  both  in  the  buckled-up  state  their  resonance  frequencies  are  close  with  a  small 
difference  due  to  fabrication  imperfections  (6.53  MHz  vs.  6.61  MHz)  as  indicated  by  the  thermal 
displacement  noise  spectrum  in  Fig.  1(b).  However,  when  one  resonator  is  displaced  from  the  up 
state  to  the  down  state,  its  mechanical  resonance  frequency  drops  to  4.05  MHz  (Fig.  1(c))  due  to 
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asymmetries  of  the  double-well  potentials.  The  intrinsic  damping  rate  of  the  resonators  y!2n  ~  2 
kHz  [13]. 
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Figure  1.  (a)  Micrograph  of  a  racetrack  cavity  with  two  llOnm  x  500nm  x  lOum  suspended  portions  as 
nanomechanical  resonators.  Insets  show  scanning  electron  micrographs  of  the  mechanical  resonators  in 
buckled  down  (left)  and  buckled  up  (right)  state.  (b,c)  thermal  noise  spectra  in  the  up-up  (b)  down-up 
state  (c).  (d)  The  measurement  setup  with  a  weak  probe  laser,  and  a  pump. 

Optomechanical  oscillators  (OMOs)  in  single  cavity,  single  resonator  systems  have  been  the 
subject  of  intense  studies  in  recent  years  [16,  17 ,  18].  It  is,  however,  expected  that  when  multiple 
oscillators  are  embedded  in  a  single  cavity  new  phenomena  will  appear  due  to  the  mutual 
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coupling  via  the  cavity  field.  Figure  2a  displays  the  evolution  of  the  RF  power  spectrum  of  the 
transmitted  probe  laser  light  (the  pump  laser  is  blue  detuned  at  A 0/k  =  0.3,  where  Ao  is  the 
detuning  between  the  laser  and  cavity  frequency,  and  k  is  the  cavity  linewidth).  At  the  lowest 
pump  powers  the  thermomechanical  motion  of  each  resonator  is  visible  as  two  lines  at  4.0  and 
6.5  MHz  respectively.  Upon  increasing  the  pump  power  in  this  first  regime,  then-  Brownian 
motion  is  amplified  as  the  backaction  reduces  the  damping.  Also,  the  optical  spring  effect  19  is 
visible  as  an  increase  in  the  resonance  frequencies.  Both  effects  are  stronger  in  the  4  MHz 
resonator  since  in  the  buckled  down  state  the  effective  refractive  index  of  the  racetrack  depends 
more  strongly  on  the  displacement  yielding  a  larger  optomechanical  coupling  gom;  From 
simulations  we  estimate  gom/27t  =  140  and  500  MHz/nm  for  up  and  down.  The  difference  in 
backaction  confirms  that  optical  backaction  is  stronger  than  photothermal  effects  as  the  latter 
would  be  the  same  on  both  resonators,  irrespective  of  their  separation  from  substrate  . 
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Figure  2.  (a)  The  evolution  of  the  RF  power  spectrum  of  the  transmitted  light  as  the  pump  power 
increases,  (b)-(e)  Cuts  through  panel  (a)  at  the  indicated  pump  power  when  both  resonators  are  in  a 
thermal  state  (regime  I,  (b)),  (c)  one  resonator  is  in  thermal  motion  while  the  other  resonator  experiences 
regenerative  oscillations  (regime  II),  (d)  the  chaotic  regime  (III)  and  (e)  the  two  resonators  are 
synchronized  (regime  IV).  The  insets  schematically  show  the  energy  of  resonators  1  (left)  and  2  (right); 
dots  correspond  to  small  thermal  motion,  and  lines  to  large  oscillations. 

When  the  pump  is  increased  beyond  -2.7dBm,  the  optomechanical  gain  fully  compensates  the 
mechanical  damping  of  the  resonator  with  the  lowest  threshold,  which  we  will  label  as  "1", 
which  starts  to  self-oscillate.  This  demarcates  the  onset  of  regime  II  which  ranges  from  -2.7  dBm 
to  0.2  dBm.  Yet,  even  though  the  oscillation  amplitude  of  resonator  1  increased  dramatically,  the 
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thermal  motion  of  resonator  2  is  undisturbed  and  is  still  clearly  visible  in  Fig.  2(a)  and  (b).  Note, 
that  the  frequency  difference  between  the  two  resonators  (2.5  MHz)  is  much  larger  than  the 
spectral  width  of  resonator  2,  so  the  cavity  occupation  oscillation  induced  by  OMO  1  cannot 
efficiently  drive  the  other  resonator.  Regime  II  thus  corresponds  to  the  single  oscillator  system 
that  we  studied  previously.  One  might  expect  that  resonator  2  simply  starts  to  oscillate  when  the 
power  is  increased  beyond  the  second  threshold  which  is  higher  due  to  its  lower  gom.  Instead, 
different  dynamics  is  encountered  (Regime  III)  where  the  power  spectrum  displays  a  large 
magnitude  over  a  wide  range  of  frequencies.  Here,  the  motion  is  chaotic  and  the  phase, 
amplitude,  and  frequency  of  the  transmitted  probe  light  change  on  short  timescales.  This  chaotic 
behavior  only  exists  for  a  limited  power  range  and  vanishes  beyond  0.7  dBm. 

Above  0.7dBm  the  two  oscillators  start  their  synchronized  motion,  as  evidenced  by  a  dramatic 
change  in  the  mechanical  displacement  spectrum.  In  Regime  IV  the  detected  photocurrent 
contains  a  strong  narrow  tone  and,  more  importantly,  the  thermomechanical  motion  of  resonator 
2  is  now  no  longer  visible.  The  single  strong  peak  in  the  RF  spectrum  indicates  that  both 
resonators  are  oscillating  at  the  same  frequency;  In  this  case  they  are  said  to  be  synchronized 
since  the  phase  difference  between  them  remains  constant.  We  have  thus  synchronized  the  two 
resonators  despite  the  extremely  large  frequency  difference:  the  second  mode,  originally  at  6.7 
MHZ,  was  almost  twice  as  fast  as  the  oscillations  of  the  first  mode  at  3.9  MHz,  indicating  the 
extremely  strong  optomechanical  interactions  and  the  tunability  of  the  double  well  potential  in 
our  system.  Note,  that  in  previous  micromechanical  synchronization  experiments  the  frequency 
difference  was  limited  to  -0.2%  [12]. 

The  spectra  in  Regime  IV  also  reveal  another  surprise:  Sidebands  emerge  around  the  strongest 
peak  (the  carrier),  indicating  a  modulation  of  the  self-sustained  oscillations.  A  close  inspection  of 
Fig.  2(a)  and  (e)  shows  that  the  spectra  contain  two  equally  spaced  sidebands  -100-500  kHz 
from  the  carrier.  Their  presence  implies  a  deteriorated  signal  phase  noise  at  that  particular 
sideband  offset  frequency  [21],  These  weak  (-20-35  dBc),  but  clearly  defined,  sidebands  are  not 
transient  phenomena  as  they  persist  during  the  entire  data-acquisition  time,  which  is  much  longer 
than  the  damping  time  of  the  resonator  y'1.  Also,  in  Regime  II  (Fig.  2(c))  with  only  a  single 
OMO  present,  sidebands  are  absent,  ruling  out  low-frequency  thermal  instabilities  [  ]  interfering 
with  optomechanical  oscillations.  Finally,  the  mixing  of  the  sidebands  in  Regime  IV  with  the 
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strong  carrier  produces  a  signal  below  500  kHz).  This  down-converted  peak  is  not  mechanical  in 
origin  but  is  due  to  the  nonlinearity  of  the  transduction  at  large  amplitudes.  Similarly,  in  Regime 
II  the  thermal  motion  of  resonator  2  is  mixed  with  that  of  OMO  1,  resulting  in  a  peak  near  2.5 
MHz. 

3.  Theoretical  analysis  of  photonic  cavity  coupled  nanomechanical  oscillators 

To  understand  the  origin  of  synchronization  and  the  slow  dynamics  in  the  cavity-coupled 
oscillators  we  theoretically  analyse  this  system;  for  details  see  [13].  Multiple  uncoupled 
oscillators  will  each  oscillate  at  their  own  frequency,  but  the  cavity  field  couples  the  oscillators 
enabling  synchronization  as  will  be  shown.  When  the  frequency  difference  between  the 
resonators  u  «  fi,  the  equation  of  motion  for  their  complex  amplitude  Uk  —  g0m,k((.uk  + 
it/c/in)exp(— iXTt) )  =  Ak  exp(i9k)  in  the  frame  rotating  at  the  average  frequency  fl  = 
i(n1  +  n2)is[i4,  23,24]: 

0i,2  =  ±i\u%2  -^f-Uli2  -  i^comX2 4>(t),  (1) 

where  uk  are  the  displacements  of  the  two  resonators  (k=  1,2).  and  <t>  (t)  describes  how  the 
photon  occupation  responds  to  a  dynamic  displacement  [24].  comk  —  hnmaxglm  k/mk[l3  are 
the  coupling  strengths  cmcl  nmax  is  the  maximum  number  of  photons  in  the  cavity.  When  multiple 
resonators  are  present,  the  cavity  only  feels  the  combined  effect  of  all  resonators.  The  same 
reasoning  as  for  a  single  OMO  shows  that  for  multiple  resonators  coupled  to  the  same  cavity,  the 
cavity  response  is  <£>(t)  =  4>(A+(t))  where  A+  is  the  magnitude  of  the  summed  complex 
amplitude  I/+(t)  =  Uk(t)  —  A+(t )  exp(i0+(t)).  A+  depends  on  the  phase  difference  between 
the  individual  oscillators,  but  not  on  the  overall  phase  9+.  The  equations  of  motion  for  the  two 
OMOs  are  thus  coupled  together  via  U+  and  the  cavity  response  function  ft>(A+).  Synchronized 
motion  of  the  two  oscillators  implies  that  they  rotate  at  the  same  frequency  fl  +  e.  Hence, 
U12(t)  —  >1,2  exp  (let)  must  be  a  solution  to  Eq.  (1):  if  no  such  solution  exists,  synchronization 
cannot  take  place.  Inserting  the  solution  into  Eq.  (1)  yields  the  equation  that  determines  the 
combined  amplitudes  of  the  limit  cycles: 
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where  com  (8com)  and  f  (5y)  are  the  average  (difference  in)  coupling  strengths  and  damping 
rates  respectively.  Solving  Eq.  (2)  can  be  done  as  illustrated  in  Fig.  3(a);  the  left  hand  side  is  a 
curve  in  the  complex  plane  parameterized  by  A+  (for  a  given  detuning  and  decay  rate),  whereas 
the  right  hand  side  depends  only  on  the  oscillators’  properties  and  is  parameterized  by  the 
unknown  frequency  e.  Intersections  of  the  two  curves  are  thus  solutions  to  Eq.  (2).  When 
inserting  the  obtained  values  for  A+  and  e  back  into  Eq.  (1)  the  individual  contributions  Y\  and  IS 
can  be  obtained  including  their  phases  relative  to  the  carrier  Y+  =  Fi+IS-  Fig.  3(b)  shows  that  for 
the  parameters  of  Fig.  3(a)  the  two  complex  amplitudes  Y\  and  IS  have  a  similar,  but  not 
identical,  magnitude  and  oscillator  1  moves  ahead  of  the  second.  Finally,  for  sufficiently 
asymmetric  oscillators,  the  two  curves  can  intersect  more  than  once,  leading  to  multistability 
[22,23],  even  in  the  unresolved  sideband  regime  where  a  single  oscillator  always  has  a  unique 
amplitude. 


Figure  3.  (a)  Complex  plane  representation  of  the  cavity  response  <t>(A+)  (black)  and  the  right  hand  side 
of  Eq.  (5)  for  two  identical  oscillators  (red)  and  two  oscillators  with  y  =  0.001  ft,  com  =  20,  v  =  0.01  ft, 
8  y  =  0.0018  ft,  and  8com  —  0  (blue).  The  curves  intersect  at  A+  =  132.0  and  e  =  0.0438  ft  (b)  Complex 
amplitudes  of  the  individual  oscillators  for  the  solution  shown  in  (a),  (c)  Colorplot  of  the  Fourier 
transform  of  the  photon  number  on  a  logarithmic  scale.  The  carrier  tone  and  its  sidebands  are  indicated  by 
red  and  white  arrows  respectively.  The  parameters  are:  Qi  =  Qi  =  6000,  k  =  526  ft,  A0  =  0.493k,  gom,i  = 

gom,2  1  - 


Equation  (2)  thus  yields  the  fixed  points  with  synchronization.  However,  to  understand  the 
dynamics  around  the  corresponding  limit  cycle,  Eq.  (1)  can  be  expanded  for  small  excursions 
and  the  eigenvalues  can  be  found.  There  are  three  independent  degrees  of  freedom:  the  two 
oscillation  amplitudes  and  the  phase  difference  between  them.  The  fourth  degree  of  freedom,  the 
overall  phase  9+,  is  not  fixed  yielding  a  zero  eigenvalue.  Depending  on  the  values  of  the  other 
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three  eigenvalues  different  types  of  behavior  are  possible:  returning  back  to  the  fixed  point 
without  oscillations  (overdamped),  returning  with  oscillations  (underdamped),  or  the  fixed  point 
is  unstable.  In  the  second  case,  when  any  of  the  two  oscillators  is  displaced,  e.g.  due  to  the  ever¬ 
present  stochastic  thermal  force  or  photon  shot  noise,  it  will  return  back  to  the  fixed  point  in  an 
oscillatory  fashion,  which  shows  up  as  sidebands  in  the  frequency  domain.  In  contrast,  for  a 
single  oscillator,  there  is  only  one  mechanical  degree  of  freedom,  the  oscillation  amplitude  A  and 
the  corresponding  eigenvalue  is  always  real.  Any  small  deviation  thus  overdampedly  returns  to 
the  limit  cycle  and  no  sidebands  are  generated  [23]. 

The  analytical  model  thus  hints  that  the  observed  sidebands  are  due  to  the  thermal  force  noise 
acting  on  the  oscillators.  To  further  analyze  the  synchronization  dynamics,  the  full  coupled 
equation  of  motion  of  the  mechanical  resonators  and  that  of  the  cavity  field  (cf.  Eqs.  (SI)  and 
(S2))  are  integrated  numerically  (see  Supplemental  Material  [13])  with  mechanical  nonlinearities 
included.  We  chose  to  simulate  a  conceptually  clear  situation  with  nearly  identical  resonators 
(5%  frequency  difference  and  identical  gom).  The  effect  of  (thermal)  force  noise  is  accounted  for 
by  kicking  the  resonators  away  from  their  steady  state  oscillations.  Figure  3(c)  shows  the 
evolution  of  the  light  field  power  spectrum  as  a  function  of  nmax  (i.e.,  the  pump  power).  Similar 
to  the  experiment  (Fig.  2),  at  low  power  two  weak  peaks  are  visible  which  also  tune  upwards 
with  increasing  power  due  to  the  optical  spring  effect.  Around  nmax  =  30  oscillations  start,  but 
now  the  two  resonators  immediately  oscillate  simultaneously  [  ].  As  expected  from  our 
analytical  theory  of  synchronized  motion,  sidebands  appear  in  the  spectrum,  but  only  when  kicks 
are  included.  As  illustrated  in  [13]  the  oscillators  are  truly  phase  locked  in  this  regime,  indicating 
full  synchronization.  When  further  increasing  nmax  the  oscillations  grow  towards  the  top  of  the 
potential  barriers  and  hence  the  frequency  goes  down.  When  the  barrier  is  crossed  at  nmax  =134 
the  detuning  suddenly  changes  dramatically  and  the  oscillations  stop.  However,  they  reappear  at 
higher  powers  and  above  the  barrier,  just  as  in  the  experiment,  the  oscillation  frequency  increases 
with  increasing  power  (see  also  Fig.  2a).  Also  note  that  the  sidebands  start  out  far  from  the 
carrier  at  the  onset  (nmax  =  335)  of  the  oscillations  in  this  regime  and  that,  just  as  in  Fig.  2a,  they 
converge  towards  the  carrier  with  increasing  power.  Finally,  the  simulations  also  reproduce 
bands  with  chaotic  behavior  with  broad  spectrums  similar  to  the  one  in  Fig.  2(d).  The 
simulations  thus  qualitatively  reproduce  most  of  the  features  observed  in  the  experiment, 
including  the  correct  tuning  behavior  of  the  resonance  frequency,  the  appearance  of  the 
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sidebands  of  the  synchronized  resonators  due  to  thermal  force  noise,  and  chaos.  It  is  expected 
that  even  better  agreement  could  be  reached  if  the  exact  shapes  of  the  mechanical  potentials  were 
known,  and  by  including  displacement-dependent  optomechanical  coupling  coefficients, 
although  we  have  to  leave  this  to  future  work. 

We  have  also  studied  the  dynamics  of  a  single  mechanical  oscillator  in  the  presence  of  an 
external  oscillator  encoded  in  the  light  field.  This  is  extremely  important  in  the  context  of 
synchronizing  a  remote  oscillator  to  an  external  clock  and  also  further  validates  our  model  for 
optomechanical  synchronization.  To  this  end,  the  pump  power  is  set  between  the  oscillation 
thresholds  of  the  first  and  second  resonator  (so  that  the  latter  does  not  play  a  role)  and  is 
modulated  at  frequency  Do  =  6.800  MHz.  When  the  modulation  index  (m)  is  zero,  resonator  1 
oscillates  freely  in  the  up-state  at  6.804  MHz  as  shown  in  the  bottom  spectrum  in  Fig.  4. 
However,  when  the  modulation  is  switched  on  the  oscillations  jump  to  Do,  synchronizing  the 
OMO  to  the  external  clock.  Interestingly,  sidebands  appear  again.  A  prominent  feature  that  the 
location  of  the  sidebands  is  not  constant:  the  offset  frequency  increases  with  m.  All  of  this  is 
reproduced  in  the  numerical  simulations  [13]  showing  that  many  of  the  phenomena  observed  in 
the  two-OMO  case  can  also  understood  in  the  conceptually  simpler  injection  locking 
experiments  [26]. 


Figure  4.  Measured  RF  power  spectral  density  of  the  detector  output  with  a  free  running  oscillator  (black) 
and  oscillations  in  the  presence  of  an  increasingly  larger  modulation  depth  of  the  pump  (red  to  dark  blue) 
for  constant  average  power.  The  curves  are  offset  for  clarity. 

4.  Prospect  and  future  plans 
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Our  technique  of  coupling  mechanical  oscillators  via  a  single  photonic  bus  creates  a  whole  new 
platform  for  nonlinear  studies.  It  will  enable  synchronization  of  large  arrays  of  individual 
optomechanical  elements  with  interesting  new  collective  phenomena  [  ]  and  allows 
synchronization  over  arbitrarily  long  distances.  Finally,  by  exploiting  the  memory  storage 
capabilities  of  the  double  well  resonators  we  envision  combining  the  information  of  the 
mechanical  bits  with  synchronization.  This  could,  for  example,  be  used  to  perform  conditional 
coupling  of  oscillators,  an  interesting  future  direction  enabled  by  our  cavity  field  coupling. 
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